Transient Zitterbewegung of charge carriers in graphene and carbon nanotubes 
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Observable effects due to trembling motion (Zitterbewegung, ZB) of charge carriers in bilayer 
graphene, monolayer graphene and carbon nanotubes are calculated. It is shown that, when the 
charge carriers are prepared in the form of gaussian wave packets, the ZB has a transient character 
with the decay time of femtoseconds in graphene and picoseconds in nanotubes. Analytical results 
for bilayer graphene allow us to investigate phenomena which accompany the trembling motion. In 
particular, it is shown that the transient character of ZB in graphene is due to the fact that wave 
subpackets related to positive and negative electron energies move in opposite directions, so their 
overlap diminishes with time. This behavior is analogous to that of the wave packets representing 
relativistic electrons in a vacuum. 



PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 
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I. INTRODUCTION 

The trembling motion (Zitterbewegung, ZB), first de- 
vised by Schroedinger for free relativistic electrons in a 
vacuum has become in the last two years subject of 
great theoretical interest as it has turned out that this 
phenomenon should occur in many situations in semi- 
conductors 0, i, i, i, i, 0, B i, m. Whenever one 
deals with two or more energy branches, an interference 
of the corresponding upper and lower energy states re- 
sults in the trembling motion even in the absence of ex- 
ternal fields. Due to a formal similarity between two 
interacting bands in a solid and the Dirac equation for 
relativistic electron in a vacuum one can use methods 
developed in the relativistic quantum mechanics for non- 
relativistic electrons in solids [ll|, Most of the ZB 
studies for semiconductors took as a starting point plane 
electron waves (see, however, Refs. [1, 0, On 
the other hand, Lock 14] in his important paper ob- 
served: 'Such a wave is not localized and it seems to be 
of a limited practicality to speak of rapid fluctuations in 
the average position of a wave of infinite extent.' Using 
the Dirac equation Lock showed that, when an electron 
is represented by a wave packet, the ZB oscillations do 
not remain undamped but become transient. In particu- 
lar, the disappearance of oscillations at sufficiently large 
times is guaranteed by the Riemann-Lebesgue theorem as 
long as the wave packet is a smoothly varying function. 
Since the ZB is by its nature not a stationary state but 
a dynamical phenomenon, it is natural to study it with 
the use of wave packets. These have become a practical 
instrument when femtosecond pulse technology emerged 
(see Ref. [l^). 

In the following we study theoretically the Zitterbewe- 
gung of mobile charge carriers in three modern materials: 
bilayer graphene, monolayer graphene and carbon nan- 



otubes. We have three objectives in mind. First, we ob- 
tain for the first time analytical results for the ZB of gaus- 
sian wave packets which allows us to study not only the 
trembling motion itself but also effects that accompany 
this phenomenon. Second, we describe for the first time 
the transient character of ZB in solids, tes ting on specific 
examples the general predictions of Ref. Third, we 

look for observable phenomena and select both systems 
and parameters which appear most promising for exper- 
iments. We first present our analytical results for bilayer 
graphene and then quote some predictions for observable 
quantities in monolayer graphene and carbon nanotubes. 



II. BILAYER GRAPHENE 

Two-dimensional Hamiltonian for bilayer graphene is 
well approximated by [l6j 



Hb = - 
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where m* — 0.054mo. The form ([T]) is valid for en- 
ergies 2 meV < £ < 100 meV in the conduction 
band. The energy spectrum is S — ±E, where E — 
h'^k'^/2m*, i.e. there is no energy gap between the 
conduction and valence bands. The position operator 
in the Heisenberg picture is a 2 x 2 matrix x{t) = 
exp{iHBt/h)xexp{—iHBt/h). We calculate 



Xii{t) 



hkH 



(2) 
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where k"^ = k"^ + k"^. The third term represents the Zit- 
terbewegung with the frequency hujz = 2h^k'^/2m*, cor- 
responding to the energy difference between the upper 
and lower energy branches for a given value of k. 

We want to calculate the ZB of a charge carrier repre- 
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sented by a two-dimensional wave packet 



V'(r,0) 



1 d 
27r 



(3) 

The packet is centered at fco — (0, k^y) and is character- 
ized by a width d. The unit vector (1,0) is a convenient 
choice, selecting the [11] component of i{t), see Ec^. 
An average of xii{t) is a two-dimensional integral which 
we calculate analytically 

x{t) = (^{r)\x{t)\^{r)) =x, + xz{t) (4) 

where Xc = {1/koy) [l — exp(— d^fepj^)] , and 



xzit) = 



«0y 



exp 



d* 



cos 



-exp(-d2fcoj^)] , 



(5) 



in which 5 — ^Jht/m* contains the time dependence. 
We enumerate the main features of ZB following from 
Eqs. Q and ([5]). First, in order to have the ZB in the 
direction x one needs an initial transverse momentum 
hkoy. Second, the ZB frequency depends only weakly on 
the packet width: ujz = {hk^y/m*){d^/{d^+S'^), while its 
amplitude is strongly dependent on the width d. Third, 
the ZB has a transient character, as it is attenuated by 
the exponential term. For small t the amplitude of xz (t) 
diminishes as exp(— F^t^) with 



Fz = 



hk. 



Oy 



(6) 



Fourth, as t (or S) increases the cosine term tends to unity 
and the first term in Eq. ^ cancels out with the sec- 
ond term, which illustrates the Riemann-Lebesgue theo- 
rem (see Ref. (3). After the oscillations disappear, the 
charge carrier is displaced by the amount Xc, which is a 
'remnant' of ZB. Fifth, for very narrow packets (d — * oo) 
the exponent in Eq. ([5]) tends to unity, the oscillatory 
term is cos{S'^kQy) and the last term vanishes. In this 
limit we recover undamped ZB oscillations. 

Next, we consider observable quantities related to the 
ZB, beginning by the current. The latter is given by 
the velocity multiplied by charge. The velocity is simply 
Vx = dxz/dt, where xz is given by Eq. (O. The calcu- 
lated current is plotted in Fig. [T}d, its oscillations are a 
direct manifestation of ZB. The Zitterbewegung is also 
accompanied by a time dependence of upper and lower 
components of the wave function. To characterize this 
evolution we define probability densities for the upper 
and lower components 



l±(3 



(7) 



where (3 = 



and the time-dependent wave 
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FIG. 1: Zitterbewegung of a charge carrier in bilayer 
graphene versus time, calculated for a gaussian wave packet 
width d — 300A and koy = 3.5 x 10*m~^: a) displacement, b) 
electric current, c) probability densities for upper and lower 
components of the wave function, d) dispersion AR{t). The 
decay time is F^^ — 40 fs, see Eq. ®. 



function is tp(r,t) — exp{—iHBt/h)ip{r,0). We have 



2tt y/n 



cos{hk^t/2m) ■ , , 
i{k+/k)^am{hkH/2m) > ' 



where k+ — k^ + iky. For t — Eq. ([5]) reduces to Eq. 
([3]). The calculated probability densities are 



P±it) 



1 1^2 

— ± — 7 exp 

2 2 s4 ^ 



d^ cos 



S^d'^E 
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where — d'^+S'^. The time dependence of P±(t) is illus- 
trated in Fig. [TJ;. Clearly, there must be P+{t) + P^{t) — 
1 at any time, but it is seen that the probability den- 
sity 'flows' back and forth between the two components. 
It is clear that the oscillating probability is directly re- 
lated to ZB. Since in bilayer graphene the upper and 
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lower components are associated with the first or second 
layer, respectively it follows that, at least for this 
system, the trembling motion represents oscillations of 
a charge carrier between the two graphene layers. For 
sufficiently long times there is P± = 1/2, so that the fi- 
nal probability is equally distributed between the two 
layers. For a very narrow packet (d — > oo) we have 
P± = (1/2)[1 ± cos^iS^k^y)], which indicates that the 
probability oscillates without attenuation. For koy = 
there is P± — (1/2)[1 ± d*/s'^], i.e. there are no oscilla- 
tions and the initial probability (1,0) simply decays into 
(1/2,1/2). 

The above phenomenon can be considered from the 
point of view of the entropy: S = —P+ lg2 P+—P- lg2 P-- 
At the beginning (the carrier is in one layer) the entropy 
is zero, at the end (when the probability is equally dis- 
tributed in two layers) the entropy is lg2 2. However, 
the entropy increases in the oscillatory fashion (see Ref. 

The transient character of ZB is accompanied by a tem- 
poral spreading of the wave packet. In fact, the question 
arises whether the damping of ZB is not simply caused 
by the spreading of the packet. To study this question 
we calculate an explicit form of the wave function given 
by integrals ([5]). The result is 



packet is not the main cause of the transient character of 
the ZB. In fact, also the spreading oscillates a little, but 
this effect is too small to be seen in Fig. [TJl. 

It is well known that the phenomenon of ZB is due to 
an interference of wave functions corresponding to pos- 
itive and negative eigen-energies of the initial Hamilto- 
nian. Looking for physical reasons behind the transient 
character of ZB described above, we decompose the total 
wave function ipir.t) into the positive [p) and negative 
(n) components ip^{r,t) and ip'^{r,t). We have 



V^"P(r,i) 



/ iyd'^koy 



■ exp 



2s4 



exp 



SW,y 

2s4 



d^ cos 



2s4 



(5^ sin 



V''°"(r-,t) 
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■ exp 



d^p" 



X exp 
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-^+d^ 

PkO 



sm 



^ PkO 
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2s4 



(11) 



where p'^ = + and pig = + [y — id'^k^yY- It is 
seen that the packet, which was gaussian at i = (see 
Eq. ([3])), is not gaussian at later times (see Discussion). 
The upper and lower components have the same decay 
time, oscillation period, etc. In order to characterize the 
spreading (or dispersion) of the packet we calculate its 
width AR{t) as a function of time 



[AR{t)r ^ {^{r,t)\r^ - {r)^mr,t)), 



(12) 



where ij}{r,t) is the above two-component wave function 
and (r) = {il}{r ,t)\f\'4){r ,t)) . The calculated width Ai? is 
plotted versus time in Fig. [TJi. It is seen that during the 
initial 80 femtoseconds the packet's width increases only 
twice compared to its initial value, while the ZB and the 
accompanying effects disappear almost completely dur- 
ing this time. We conclude that the spreading of the 



\m) 



-iHt/h 



\m) 



e-^'^'^^imom + e'^*/''(n|V(0))|n), (13) 



where \p) and \n) are the eigen- functions of the Hamilto- 
nian ([1]) in fc space corresponding to positive and negative 
energies, respectively. Further 



{k\p) = ^ ( ,2^/1.2 ]S{k-k ), 



(kin) 



V2 \ -kl/k 



2 /,^2 ]5{k-k 



(14) 
(15) 



After some manipulations we finally obtain 



(10) 



The function 'ip"{r, t) is given by the identical expression 
with the changed signs in front of E and fc^/fc^ terms. 
There is ■0(^,0 = ip^ {r , t) + {r , t) and (■0"|V'^) = 0. 

Now we calculate the average values of x and y using 
the positive and negative components in the above sense. 
We have 



x{t)^ / (V'" + V^)^a;(V'" + V'^)d' 



(17) 



so that we deal with four integrals. A direct calculation 

gives 



\^xdh 



(18) 



i;''''xi;Pd'^r + i;P''xij"'d'^r = xz{t), (19) 



where Xc and xz{t) have been defined in Eq. ([?]). Thus 
the integrals involving only the positive and only the neg- 
ative components give the constant shift due to ZB, while 
the mixed terms lead to the ZB oscillations. All terms 
together reconstruct the result (g]). 

Next we calculate the average value y. There is no 
symmetry between x and y because the wave packet is 
centered around fc^; = and ky = k^y . The average value 
y is again given by four integrals. However, now the 
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mixed terms vanish since they contain odd integrands of 
kx , while the integrals involving the positive and negative 
components alone give 



2m* 
2m* ' 



(20) 
(21) 



This means that the 'positive' and 'negative' subpackets 
move in the opposite directions with the same velocity 
V = hkQyt/2m*. The relative velocity is w'"'^' = Hkoyt/m*. 
Each of these packets has the initial width d and it 
(slowly) spreads in time. After the time T"^^ — d/v^*^^ 
the distance between the two packets equals d, so the 
integrals (jlSp are small, resulting in the diminishing Zit- 
terbewegung amplitude. This reasoning gives the decay 
constant Tz = hkQyt/m*d, which is exactly what we de- 
termined above from the analytical results (see Eq. ([6])). 
Thus, the transient character of the ZB oscillations is 
due to the increasing spacial separation of the subpackets 
corresponding to the positive and negative energy states. 
This confirms our previous conclusion that it is not the 
packet's slow spreading that is responsible for the atten- 
uation (see Discussion). However, as we show below, also 
spreading may possibly play this role in some cases. 

To conclude our analytical discussion of the ZB in bi- 
layer graphene we consider an interesting property of the 
velocity squared. If Vx = dH/dpx and Vy — dH/dpy are 
calculated directly from the Hamiltonian ([1]), then it is 
easy to show that = H'^k'^/m'^ and Vy — h^k'^/m'^, so 
that = v'^+Vy = 2h'^k'^/m^ does not depend on time. 
In the Heisenberg picture we split the velocity compo- 
nents into 'classical' and ZB parts 



_ hky / k"^ s\n{hk'^t/m) 



ik"^ cos{hk'^t/m) 



-ik"^ cos{hk'^t/m) 
-k^ sm{hk^t/m) 







kl 




(22) 



and similarly for Vy{t). Noting that {wf (^), w^(i)} = 
we have Vx{tf = v^{tf + v%{tf . Using Eq. ^ 
we show that each of these terms is time independent: 
^x{t)'^ = h^ky/rri^ and v%{tY = h^kl/rri^, and simi- 



larly for Vy{tY and Vy(t)'^. 



Thus, the velocity squared of 



the ZB component v^{tY=v^{tY + {ty =h^ k^ / is 
equal to that of the 'classical' component v'^ {tY =v'!^{tY + 
v!y{tY=h^k^/m^. 

III. MONOLAYER GRAPHENE 

Now we turn to monolayer graphene. The two- 
dimensional band Hamiltonian describing its band struc- 
ture is M M,M,M,M 



Hm = u 





Px + iPy 



Px - Wy 





(23) 



where m « 1 x lO^cm/s. The resulting energy disper- 
sion is linear in momentum: £ = ±uhk, where k = 

kl + fc^ . The quantum velocity in the Schroedinger 

picture is Vi = OHm / dpi, it does not commute with the 
Hamiltonian (j23p . In the Heisenberg picture we have 
v{t) — exp{iHMt/h)vexp{—iHMt/h). Using Eq. 
we calculate 



■ sm{2ukt). 



(24) 



The above equation describes the trembling motion with 
the frequency ljz = 2uk, determined by the energy differ- 
ence between the upper and lower energy branches for a 
given value of k. As before, the ZB in the direction x oc- 
curs only if there is a non- vanishing momentum hky. We 
calculate an average velocity (or current) taken over the 
wave packet given by Eq. ([3]). The averaging procedure 
amounts to a double integral. The latter is not analytical 
and we compute it numerically. The results for the cur- 
rent jx = evx are plotted in Fig. [2]for koy = 1.2 x 10®m~^ 
and different realistic packet widths d (see Ref. [1^). It 
is seen that the ZB frequency does not depend on d and 
is nearly equal to ojz given above for the plane wave. On 
the other hand, the amplitude of ZB does depend on d 
and we deal with decay times of the order of femtosec- 
onds. For small d there are almost no oscillations, for 
very large d the ZB oscillations are undamped. These 
conclusions agree with our analytical results for bilayer 
graphene. The behavior of ZB depends quite critically 
on the values of k^y and d, which is reminiscent of the 
damped harmonic oscillator. 



IV. CARBON NANOTUBES 

Finally, we consider monolayer graphene sheets rolled 
into single semiconducting carbon nanotubes (CNT). 
The band Hamiltonian in this case is similar to Eq. P5)) 
except that, because of the periodic boundary conditions, 
the momentum px is quantized and takes discrete values 
hkx = hknu, where km, = (27r/L)(rt — 1^/3), n ~ 0, ±1, . . ., 
z/ = ±1, and L is the length of circumference of CNT 
(2^ [27t. As a result, the free electron motion can occur 
only in the direction ?/, parallel to the tube axis. The ge- 
ometry of CNT has two important consequences. First, 
for V — ±1 there always exists a non- vanishing value of 
the quantized momentum hknu- Second, for each value 
of knu there exists k^n.-u = ^km, resulting in the same 
subband energy £ = ±E, where 



E = huxlk'^^^^ 



k^. 



(25) 



The time dependent velocity Vy{t) and the displacement 
y{t) can be calculated for the plane electron wave in 
the usual way and they exhibit the ZB oscillations (see 
'"I). For small momenta ky the ZB frequency is 
2huknu- The ZB amplitude is 



Ref. 

hioz 



Eg, where Eg 
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2 4 6 8 t [fs] 12 

FIG. 2: Oscillatory electric current caused by the ZB in 
monolayer graphene versus time, calculated for a gaussian 
wave packet with koy = 1.2 x 10®m~^ and various packet 
widths d. 



~ ^/knu- However, we are again interested in the dis- 
placement y{t) of a charge carrier represented by a one- 
dimensional wave packet analogous to that described in 
Eq. (HD 

^(y) = J dkye-y'^ie^'^^y J ) . (26) 

The average displacement is y{t) = yz{t) — Vsh where 
, , h^du^hau dky f2Et\ .2 ,,2 

and ysh = l/20r(i sgn(6)[l - exp(62), where 

b = knud and $(a;) is the error function. The ZB oscilla- 
tions of y{t) are plotted in Fig. [3] for n — v — ±1 
and L = 200 A. It can be seen that, after the tran- 
sient ZB oscillations disappear, there remains a shift ysh- 
Thus the ZB separates spatially the charge carriers that 
are degenerate in energy but characterized by n, v and 
— n, —V quantum numbers. The current is proportional 
to Vy = dy/dt, so that the currents related to = 1 
and u ~ ~\ cancel each other. To have a non- vanishing 
current one needs to break the above symmetry, which 
can be achieved by applying an external magnetic field 
parallel to the tube axis [J] . 

Above we considered the situation in which a non- 
vanishing value of transverse momentum hkx is 'built 
in' by the tube topology. However, it is also possible 
to prepare a wave packet with an initial non-vanishing 
momentum fcoy Using the method presented above for 
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1 10 100 t [fs] 1000 

FIG. 3: Zitterbewegung of two charge carriers in the ground 
subband of a single carbon nanotube of L = 200 A versus 
time (logarithmic scale), calculated for gaussian wave packets 
of two different widths d. After the ZB disappears a constant 
shift remains. The two carriers are described by different 
quantum numbers v. At higher times the amplitude of ZB 
oscillations decays as t~^^'^ (see text). 

bilayer graphene (see Eq. (fT5|) ) we can decompose the 
total wave packet (|26p into the positive and negative sub- 
packets with the result 

(28) 

The function ■0" (y, t) is given by a similar expression with 
the changed signs in front of E and (fc+ /k) terms. Here 

we use the notation k = y^fc^v^^TA?^ and fc+ — knu + iky. 

Now the oscillating part of y is, as before 

j r^y^dy + J r^yrdy = yz{t, koy). (29) 

For koy — the above yzit,koy) reduces to yzit) given 
by Eq. ([TT)) . The average contributions of positive (or 
negative) terms alone are 

/ - ±ut[ fc'^l^fel^. , (30) 

2 'i^kl^ + kl 

where F{ky) = d^/'^ /{2Tr^/^) exp{-^d^{ky - k^yf) is the 
packet function. The sum of the first terms for tljP{y,t) 
and ^^{y,t) in Eq. ([50)1 gives ijc, as before. For koy = 
the second term vanishes which physically means that 
the relative velocity of the two subpackets is zero, so that 
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they stay together in time. It is for this reason that the 
decay of ZB is slow (see Fig. If fcgy 7^ 0, the second 
term in Eq. ([30| does not vanish, the two subpackets run 
away from each other, their overlap diminishes and the 
ZB disappears much more quickly. 

The question remains: what is the physical reason for 
the (slow) damping of the ZB electron shown in Fig. [3l 
if the subpackets stay together? (As we mentioned in 
the Introduction, the mathematical expression for the 
damping phenomenon is the Riemann-Lebesgue theo- 
rem.) Trying to answer this question we calculated the 
spreading of the wave subpackets in time. For the 
initial width Ay 90 A the subpackets reach the width 
Ay 2600 A after the time of 1000 fs (see Fig. ^. 
Thus, we would be tempted to say that it is the spread- 
ing of the packets that is responsible for the attenuation 
of ZB. However, it should be noted that, while at higher 
times the packet's dispersion is linear in time (see Ref. 

and Fig. [3li), the amplitude of ZB oscillations de- 
cays as t~^/^. A similar slow damping of ZB occurs for 
one-dimensional relativistic electrons in a vacuum if the 
average momentum of the subpackets is zero (see Discus- 
sion). 



V. DISCUSSION 

It is of interest that the ZB phenomena similar to those 
described above occur also for wave packets representing 
free relativistic electrons in a vacuum governed by the 
Dirac equation. This confirms again the strong similar- 
ity of the two-band models for non-relativistic electrons 
in solids to the description of free relativistic electrons in 
a vacuum, see Refs. [1, 0, 0, [lil, US]- In contrast to bi- 
layer graphene, the kinematics of the one-dimensional rel- 
ativistic wave packets may not be described analytically, 
so the solutions were computed numerically and visual- 
ized graphically by Thaller ^2^. It was shown that: 1) An 
initial relativistic gaussian wave packet after spreading is 
not gaussian any more. This is analogous to our Eqs. 
pO)) and (|lip . 2) If an average momentum of the initial 
positive and negative subpackets is zero, the overlap of 
the two subpackets remains almost constant in time and 
the resulting ZB decays quite slowly. This corresponds to 
our considerations of CNT with koy = 0, see Fig. [3l (It 
is to be reminded that the two overlapping subpackets 
are orthogonal to each other.) 3) If the initial average 
momentum of both subpackets is nonzero, the two sub- 
packets quickly run away from each other and the ZB 
falls quickly since it is sustained only when the subpack- 
ets have some overlap in the position space. This corre- 
sponds to our considerations of bilayer graphene, see Fig. 

HI 

The transient ZB of free relativistic wave packets in 
a vacuum was also studied numerically by Braun et al. 
[2^. It was shown that, for example, the decay times 
of a typical wave packet having the width Aa; = Ac and 
the initial wave vector fcox = 1.37 a.u. is 2.4 x 10~^ fs. 



This should be compared with our predicted decay times 
of = 40 fs for bilayer graphene. It turns out once 
again that solids are much more promising media for an 
observation of Zitterbewegung than a vacuum. 

The Zitterbewegung phenomenon described above 
should not be confused with the Bloch oscillations of 
charge carriers in superlattices, although the latter occur 
at picosecond frequencies and have comparable picosec- 
onds decay times (see e.g. Refs. [13, Hi], H^])- How- 
ever, the Bloch oscillations are basically a one-band phe- 
nomenon, they have been realized in superlattices (al- 
though this is in principle not the condition sine qua non) 
and, most importantly, they require an external electric 
field driving electrons all the way to the Brillouin zone 
boundary. On the other hand, the ZB needs at least two 
bands and it is a no-field phenomenon. On the other 
hand, narrow-gap superlattices could provide a suitable 
medium of its observation. 

In view of our results it is clear that, in order to observe 
the transient Zitterbewegung, it is necessary to prepare 
simultaneously a sufficient number of charge carriers in 
the form of wave packets. If one wants to detect the 
current, the trembling motion of all carriers must have 
the same phase. On the other hand, if one wants to see 
only the remnant displacement, the phase coherence is 
not necessary. As we said above, the ZB frequency is to 
a good accuracy given by the corresponding energy differ- 
ence between the upper and lower energy branches while 
the amplitude depends strongly on packet's width. For 
the two graphene materials considered above one needs 
an initial momentum in one direction to have the ZB 
along the transverse direction (see also Ref. Q). For 
nanotubes the initial momentum is automatically there 
due to the circular boundary conditions. The oscillatory 
motion between two graphene layers, as illustrated in Fig. 
[1]:, appears to be a promising phenomenon for observa- 
tion. As far as the detection is concerned, one needs 
sensitive current meters or scanning probe microscopy, 
both working at infrared frequencies and femtosecond to 
picosecond decay times (see Refs. [33l.[3^). 



VI. SUMMARY 

In summary, using the two-band structure of bilayer 
graphene, monolayer graphene and carbon nanotubes we 
show that charge carriers in these materials, localized in 
the form of gaussian wave packets, exhibit the transient 
Zitterbewegung with the decay times of femtoseconds in 
graphene and picoseconds in nanotubes. Observable dy- 
namical ZB effects, most notably the electric current, are 
described. It is demonstrated that, after the trembling 
motion disappears, there remains its 'trace' in the form of 
a persistent charge displacement. It is emphasized that 
the described ZB in solids is in close analogy to that of 
the relativistic electron in a vacuum. 
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